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The cubic kilometer IceCube neutrino telescope now operating at the South Pole in a near 
complete configuration observes the neutrino sky with an unprecedented sensitivity to galactic 
and extra- galactic cosmic ray accelerators. Within the multi-messenger framework, IceCube 
offers unique capabilities to correlate and contrast the neutrino sky with the gamma-ray sky 
and ultra high energy cosmic rays and complements other indirect and direct dark matter 
search programs. We review here the status of the experiment and recent selected results. A 
discussion of the implications of the observations will be followed by the prospects of future 
developments, substantially extending the reach of the observatory at extremely high energies, 
in the GZK region and at low energies enhancing capabilities to probe dark matter and cosmic 
ray sources in the southern sky. 



I. Introduction 

The IceCube construction signals the emergence of a new class of gigantic detec- 
tors dedicated to the observation of the high energy (HE) neutrino sky. IceCube's 
main goals are the unambiguous identification of the first galactic and extra-galactic 
cosmic ray (CR) accelerators with the detection of HE neutrinos from point sources 
and the divulgation of the nature of dark matter (DM) through the observation of 
a secondary neutrino flux from annihilating DM in our galaxy. A multi-messenger 
approach applies particularly well to these research topics: the knowledge which 
can be gained by combining and contrasting measurements by means of various 
astroparticles strongly enhances the physics return for the community. 

We introduce the IceCube neutrino telescope in section [Hi and present selected 



results revealing the multi-messenger approach potential in section III : Searches 
for DM annihilation in the galactic halo and for neutrinos from the X-ray bi- 
nary LS 1+61° 303 based on time-dependent neutrino flux predictions in sec- 



tions |III A| and |IIIB a search for a correlation of HE neutrinos with ultra high 



energy (UHE) CR and the prospects for UHE neutrino detection in sections IIIC 



and III D| . IceCube as a probe for CR anisotropics in the Southern sky is presented 



in section IIII El 



II. The IceCube detector 

The current layout of the IceCube neutrino observatory consists of an in-ice ar- 
ray of 79 strings deployed at a maximal depth of 2.45 km under the South Pole ice 
cap instrumenting a km 3 of crystal clear ice and a surface air shower array IceTop 
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instrumenting a square km. Since the deployment began in 2004, the detector has 
been operating in the various configurations, referred to as IC-xx, where xx is the 
number of strings in operation, reflecting the construction status. The observatory 
is scheduled for completion in December 2010, after a very successful deployment 
over the past five years [Tj. The design is meeting performance expectations [2]. The 
in-ice array will eventually be equipped with 86 strings, a nominal string spacing 
of 125 m, with a maximal number of 60 digital optical modules (a photomultiplier 
and electronics for signal digitization, time-stamping and communication with the 
ground surface) [31 H] . In the center of the in-ice array, 6 strings more densely instru- 
mented form together with the neighbouring strings a dense inner core, enhancing 
the IceCube detection capabilities toward lower energies and potentially enabling 
IceCube to explore the Southern neutrino sky. 

In its current state, IC-79 is taking data at a rate of approximately 2 kHz, dom- 
inated by secondary muons originating from cosmic air showers in the atmosphere. 
The current integrated exposure is equivalent to about two years with the com- 
pleted IceCube detector. Moreover, with the now operational deep core, the data 
are enriched with lower energy neutrinos. After the application of filtering and re- 
construction procedures to the events [5] , the major experimental muon background 
can be largely reduced in order to keep a sample of events dominated by neutrino- 
induced muon events. Astrophysical searches look for an excess over atmospheric 
neutrino background expectations. 

The IceTop air shower array is dedicated to composition studies of the CR spec- 
trum around the knee and above [BJ. IceTop enables the detection of air showers 
at energies above ~1 PeV. As a part of IceCube, it also plays an important role 
of facilitating background rejection and calibration of events triggered by the in-ice 
array [7]. 

On site, prospecting activities for the detection of UHE neutrinos are taking place 
as well, relying on the alternate radio Cherenkov and acoustic signatures accompa- 



nying UHE neutrino interactions, see Set. HID 



III. Multi-messenger Astroparticle Physics with IceCube 
A. Searches for supersymmetric dark matter 

Under certain assumptions, the minimal supersymmetric standard model 
(MSSM) framework provides a stable weakly interacting DM candidate: the neu- 
tralino, a self-annihilating thermal relic of the early universe. Its mass, bounded 
from below by accelerator constraints and above by theory, lies between 46 GeV [8] 
up to a few TeV |9J. Secondary particles, including is, are emitted at a higher 
rate from regions of greater DM density where gravitationally trapped neutralinos 
annihilate pairwise [10]. 

The galactic halo or compact objects such as the Sun seem to be promising regions 
of such enhanced DM densities and thus for conducting dedicated searches for these 
signatures. These analyses all search in common for an excess from the directions 
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of these enhanced DM density regions. In the absence of any excess over the known 
atmospheric neutrino background, upper limits on the neutrino-induced muon flux 
from DM annihilations are obtained. The deeper connection to the physics arises 
with the conversion of the neutrino-induced muon flux limits into cross section upper 
limits [llj: self-annihilation cross sections, velocity averaged (in halo analyses), cr^v, 
and spin-dependent scattering cross sections ctsd in search for signatures from the 
self-annihilating solar neutralinos (assuming equilibrium between capture and anni- 
hilation rate in the Sun). Assuming neutralinos constitute a sizable fraction of the 
galactic DM density, these analyses have a significant potential to exclude regions 
from the MSSM parameter space which would otherwise remain unconstrained by 
direct search experiments [12J and by 7-ray and CR experiments performing indirect 
searches similar to IceCube [13] . 



90% C.L. Upper Limit - Einasto %X~ >TT 




FIG. 1: Left: limits on a aV w.r.t. m x for various self-annihilation channels and both halo 
analyses, including SK results for comparison [32 [53]. Right: limits of the halo analyses 
and favored region by PAMELA and Fermi for the rr self-annihilation channel. 

Novel results from searches for a DM signature from the halo are presented here 
as they seem promising; indeed powerful as well as competitive and complementary 
with other indirect searches. Two different searches were performed, respectively 
with IC-22 and IC-40 [Hj. The first search looked for a differential excess from 
the northern sky selecting two broad regions, one pointing toward the Galactic 
Center (but not including it, as it is located in the Southern hemisphere) and a 
corresponding mirror region (defined by a shift of 180°in right ascension). This 
analysis is referred to as the outer Galaxy analysis. The second and more ambitious 
is a search for an excess from the Galactic Center (GC); a difficult undertaking due 
to the heavy atmospheric muon background contamination as the GC is above the 
horizon (i.e. in the Southern hemisphere). The search is made possible for the first 
time with IC-40 thanks to the conjunction of a large detection volume enabling to 
veto the background efficiently and the fact that the signal, according to specific 
models profiling the DM density in the Milky Way [15 1, is strongly increasing in 
this direction (the neutralino annihilation rate scales with the square of the DM 
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density [IS]). This analysis is however slightly more dependent on the assumed DM 
profile, uncertain near the GC. 

Results of these IceCube's searches are presented in Fig. [T] From a multi- 
messenger perspective they are remarkable: combining the measurements from the 
Fermi and PAMELA experiments and in a leptophilic (xx ~^ TT ) neutralino self- 
annihilation scenario |17| . the favored region is excluded by IceCube's GC analysis 
(considering the Einasto profile). This is a compelling illustration of the great con- 
trasting potential of observations, which is achieved by combining observational 
constraints issued by means of the various messengers. Note however, that the 
direct superposition on a single plot of results, which rely on the validity of the 
assumed DM profile, must be interpreted with care: the neutrino flux exclusion 
limit comes from the integrated self-annihilating DM rate along the line of sight, 
while the anomalous measurements of the positron fraction by PAMELA and of the 
electronic component by Fermi signs local annihilation, as electrons and positrons 
cool off quickly. A more conservative interpretation of the anomalous measurements 
calls for local CR accelerators [T8H2T] . 

B. Search for neutrinos from the X-ray binary LS 1+61° 303 

While the resolution of the centennial mystery of CR origin constitutes one of 
IceCube main's goals as it would deepen our understanding of the astrophysical 
sources which accelerate them, hadronic accelerators have not yet been revealed. 
Tight constraints have instead been issued with an average upper flux limit from a 
random excess in the northern sky which lies at E 2 dQ/dE ~ 0.5 • 10 -11 TeV cm -2 
s~ 1 (IC-40) for Fermi- accelerated CR in steady point source of neutrinos |24| . This 
limit can be substantially improved however at the expense of an increased model 
dependency of the search, for instance by looking at an excess from predefined 
sources or source class (source stacking) in a catalog [2U [25] . This limit can be 
further overcome by accounting for time-dependence of the neutrino emission in the 
test statistics with an assumed neutrino flux enhancement related to the variability 
of the multi- wavelength spectrum [26]: 

1. If the source is not periodic and exhibits random MWL variability, e.g. blazar, 
analyses assuming a positive correlation between specific photon bands and the 
neutrino flux were conducted [23]. These searches are strongly model-dependent. 

2. If the source is periodic, e.g. the gamma-ray loud X-ray binary 
LS 1+61° 303 considered below, we have considered two approaches: either as- 
suming an equally periodic enhancement of the neutrino emission, but unknown 
normalization and phase as described in |27| or, alternatively, applying a new anal- 
ysis methodology, based on specific models relating multi- wavelength spectra to the 
expected spectral neutrino "light curves" [28]. While this second approach suffers 
a reduced statistical penalty (the price to pay for the fit of the parameters) and 
may thus boost the IceCube's discovery potential, it makes the search more strongly 
model-dependent. 

We explore here this new analysis methodology and apply it to the fascinating 
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gamma-ray loud X-ray binary LS 1+61° 303, a periodic neutrino source candidate 
using the hadronic model described in |29l [30] among others |31H35| and recently 
subject to intense multi- wavelength observation campaigns [36-40J. Together with 
LS 5039 and PSR B1259-63, located in the Southern hemisphere and therefore out 
of the field of view of IceCube, the TeV emitter LS 1+61° 303 exhibits a strong com- 
ponent in the GeV band, variability over two timescales, the orbital period with 26.5 
days and the superorbital period with 4.6 yrs, also a notable source of outburst [4T] . 
The nature of the compact object is not yet known |34j and the related inclination of 
the system orbital plane w.r.t. the line of sight not well constrained |42j . The mas- 
sive companion is a Be star with a surrounding decretion disk periodically disrupted 
over the superorbital timescale, suggesting a possible modulation of the neutrino 
emission. Near the apastron, the source exhibits enhanced VHE, X-ray and radio 
activity |36[ [371 14^1 144] anti-correlated with the GeV emission [45] . 



LS I +61 303: Ephemeils of the Radio Outburst Modulation and AM AND A/Ice Cube event candidates 
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FIG. 2: Left: selected neutrino candidate events with AMANDA, IC-22 and IC-40 and the 
superobital radio modulation within the emission window of a model (see text). Right: 
phase and energy estimate of the selected events in IC-22. 

Following an analysis optimizing the discovery potential, candidate neutrino 
events were selected from IC-22 and IC-40. An event subsample from the publicly 
released AMANDA events |46j was selected using a similar analysis. Fig. [2j left, 
shows the final event sample recorded by AMANDA, IC-22 and IC-40 in the last 10 
years together with the radio modulation. On the right, the phase of the IC-22 se- 
lected events is presented. The a-posteriori p- value of about 1% is obtained from the 
HE AMANDA events. The p-value of the IC-22 analysis is approximately 2%. The 
p- value extracted from the blind IC-40 search is not significant. Refined models for 
neutrino flux predictions including a mechanism to explain the superorbital phase 
(such as a Be star disk growth followed by a disruption) and the radio and X-ray 
variabilities on this time scale could be of interest for a long term data analyses 
of the neutrinos from the direction of LS 1+61° 303. In conclusion, this approach 
set more stringent constraints on neutrino emission models and may facilitate the 
discovery and characterization of a neutrino signal from galactic and extra-galactic 
point sources. Upon discovery, not only a first incontrovertible demonstration of 
hadronic acceleration in the environment of specific astrophysical systems would be 
provided, but the underlying source acceleration mechanisms and morphology would 
be constrained. 
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C. Ultra high energy cosmic rays and neutrinos 

The Pierre Auger Observatory (PAO) and HiRes experiments have in the past 
decade reported 35 events above 57 EeV, pointing back to the northern hemi- 
sphere [27]. At these high energies, UHE CRs likely travel only weakly deflected 
over long extra-galactic distances, thus approximately point back to the cosmic 
accelerator from which they originate. Therefore, associated neutrinos produced ei- 
ther locally or along the propagation are closely aligned, motivating a search for a 
directional correlation between neutrinos and UHE CRs. Moreoever the strong cor- 
relation of the UHE CR arrival directions with nearby active galactic nuclei which 
was earlier reported by PAO makes this analysis exciting. With the release of an 
updated UHE CR catalog [38], the association now has a reduced significance, but 
there are still indications for the possible existence of such a correlation. 

An analysis was conducted with a subsample of the IC-22 data, enriched with 
HE neutrinos (with a mean energy 100 TeV for an E~ 2 simulated spectrum): the 
neutrino candidate events were selected on the basis of their angular distance \l/ < 3 ° 
to one of the 35 reported UHE CR events. 60 events were found, while 43.7 events 
were expected ("off-source" estimation). The probability of such an excess is 0.98% 
in the background-only hypothesis. A similar analysis was repeated with the IC- 
40 data, including the events from the updated PAO catalog, bringing the number 
of UHE CR events to 82: the previous excess is partially washed out (298 events 
observed, 274 expected). 

D. Alternative GZK neutrino detection technique 

The IceCube potential for the observation of GZK |49] neutrinos strongly de- 
pends on the GZK flux normalization, which depends on source evolution, injection 
spectrum and CR composition [50] . The characterization of the GZK neutrino flux 
would enable the partial recovery of the degraded information carried by UHE CR, 
which would permit the delineation of cosmological source evolution scenarios from 
source injection spectrum characteristics. This in turn would help determine the 
nature of the most powerful CR accelerators in the universe. 

Currently, the situation concerning the GZK neutrino flux normalization is un- 
certain [51} 152] : while the observed correlation of UHE CR sources with the AGN 
distribution by AUGER [47J hints at a light composition (and in this case current 
neutrino flux limits lie close to the upper flux predictions [53J), dedicated AUGER 
composition studies favor a composition becoming heavier at UHE [53]. It is there- 
fore necessary to guarantee an observation of GZK neutrinos (and perhaps to char- 
acterize its flux) to build detectors with substantially larger effective volumes than 
IceCube. However, these detectors should not rely on the optical signature, as 
this detection technique has practical limitations linked to a rather short attenua- 
tion length for Cherenkov light propagation in the ice [55]. Alternative detection 
techniques are called for [56-61J: A rather promising technique is based on the ob- 
servation of the radio Cherenkov signature of the interaction shower [H2] • Over the 
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past decade the ice has been well characterized with the RICE detector and the 
attenuation length of the radio signal was demonstrated to be much longer than for 
the optical signal [63], correspondingly enabling a sparser instrumentation of the 
detection volume. The Askaryan Radio Array (ARA) [M] project is now approved 
for a 3-year preparatory phase, aiming at demonstrating the feasibility of the full 
scale detector from technological and environmental aspects. ARA would eventually 
be instrumenting 80 km 2 and be providing a 3-year sensitivity sufficient for the de- 
tection of GZK neutrinos induced by photo-disintegration in pessimistic predictions 
of pure iron composition. An alternative technique relies on detection of the acous- 
tic signature of the GZK neutrinos. Its applicability is studied by the South Pole 
Acoustic Test Setup (SPATS) [65] . It was recently found that the attenuation length 
is not as long as anticipated [661 EZ] , challenging the design of the acoustic detection 
devices [68] . The experimental effort will continue with the deployment of additional 
devices this winter at the South Pole aiming to measure the absolute level of envi- 
ronmental noise, which may, if too high, impede the applicability of the technique in 
this location. The technique, yet to mature from technological and environmental 
characterization standpoint, could eventually be used in conjunction with the radio 
detection technique: the detection of a hybrid event, i.e. detected with two or more 
distinct signatures (optical, radio or acoustic), would unambiguously demonstrate 
the interaction of a UHE neutrino. 



E. Cosmic ray anisotropies 

Diverting the multi-messenger approach into multi-messenger detection capabil- 
ities of the IceCube's in-ice component, we have searched and found first evidence 
for southern sky anisotropies in the arrival direction of CRs. Analyses of the IC-22 
and IC-40 data are described in detail in [691 [70] . The sample of IC-22 down-going 
muon data reported at this conference shows an anisotropy. The maximal peak to 
peak amplitude of the projected skymap in right ascension is about 0.15%. More- 
over, when considering two event subsamples with increasing median energies, from 
12.6 to 126 TeV, the relative large scale anisotropy reduces significantly. The com- 
bined significance skymap of TeV MILAGRO and mufti- TeV IC-40 events is shown 
Fig. [3] suggest that these medium scale anisotropies could be part of a larger scale 
anisotropy. 

The origin of these anisotropies, reported by several surface and underground 
experiments located in the northern hemisphere |71|, [72] in the past two decades, 
is not yet well understood: anisotropies are expected due to the Compton-Getting 
effect [73], but have not yet been confirmed. On the contrary, the anisotropy seems 
to be out of phase with this effect. The tail- in excess combined with the energy 
dependence of the observed anisotropy may suggest a local origin, explained by 
heliospheric effects [73] , or other local features of the magnetic field over larger scales, 
which would enable the particles to stream from nearby cosmic accelerators. In this 
case, it would as well establish a connection with the electronic component anomalies 
discussed in section Ull A[ It is clear from the above that the better characterization 
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FIG. 3: Preliminary sky map significance combining data from IC-40 and MILAGRO. 

of the anisotropies, at various angular scales, could be quite revealing of the involved 
underlying physics. 



IV. Conclusion 

We have presented recent results at the interface of the three pillars of astropar- 
ticle physics: particle physics, astrophysics and cosmology, with a clear focus on 
the multi-messenger approach, which may help delineate the nature of the cosmic 
accelerators or unveil the nature of the DM. 

IceCube has recently realized its potential for the study of transient events: the 
IC-40 neutrino flux limits on Gamma- Ray Bursts (see e.g. [TDJ [75] ) is now excluding 
the Waxman-Bahcall upper bound [76] and in the future, fully exploiting the deep 
core in the multi-messenger context, IceCube will certainly be revealed as a powerful 
tool for issuing strong constraints on the nature of the DM and for the observation 
of MeV neutrinos from galactic supernova [77J [78] • 



Acknowledgments 

The author, supported by the Swiss National Research Foundation, grant No 
PP002 114800, warmly thanks the organizer for a great conference. 



[1] R. Abbasi et al, NIM A 601 (2009) 294; M. Ackermann et al, ApJ 675 (2008) 1014. 

[2] A. Achterberg et al. [IceCube coll.], ApP 26 (2006) 155. 

[3] R. Abbasi et al. [IceCube coll.], NIM A 601 (2009) 294. 

[4] R. Abbasi et al. [The IceCube coll.], NIM A 618 (2010) 139. 

[5] J. Ahrens et al. [AMANDA coll.], NIM A 524 (2004) 169; G. Japaridze, M. Ribordy, [astro- 
ph/0506136]. 

[6] J. Ahrens et al. [AMANDA & SPASE coll.], ApP 21 (2004) 565. 

[7] J. Ahrens et al. [SPASE & AMANDA coll.], NIM A 522 (2004) 347. 

[8] J. Abdallah et al. Eur. Phys. J. C 31 (2003) 421. 

[9] K. Griest, M. Kamionkowski, Phys. Rept. 333 (2000) 167. 

[10] I. F. M. Albuquerque and C. Perez de los Heros, PRD 81 (2010) 063510. 



9 



[11] G. Wikstrom, J. Edsjo, JCAP 0904 (2009) 009. 

[12] F. Halzen and D. Hooper, New J. Phys. 11 (2009) 105019. 

[13] R. Abbasi et al [IccCubc coll.], PRL 102 (2009) 201302. 

[14] C. Rott for the IccCubc coll., arXiv:0912.5183 [astro-ph.HE]. 

[15] B. Moore et al, MNRAS 310 (1999) 1147; J. F. Navarro, C. S. Frcnk and S. D. M. White, ApJ bf 

462 (1996) 563; A. V. Kravtsov et al, ApJ 502 (1998) 48. 

[16] H. Yuksel, S. Horiuchi, J. F. Beacom et al, PRD 76 (2007) 123506. 

[17] P. Meade et al, Nucl. Phys. B 831 (2010) 178. 

[18] D. Hooper, P. Blasi, P. D. Serpico, JCAP 0901 (2009) 025. 

[19] H. -B. Hu, Q. Yuan, B. Wang et al, ApJ 700 (2009) L170. 

[20] M. Ackermann et al [Fermi LAT coll.], arXiv:1008.5119 [astro-ph.HE]. 

[21] H. Yuksel, M. D. Kistlcr and T. Stancv, PRL 103 (2009) 051101. 

[22] S. Dcsai et al [SK coll.], PRD 70 (2004) 083523. 

[23] J. Hisano, K. Nakayama, M. J. S. Yang, Phys. Lett. B 678 (2009) 101. 

[24] R. Abbasi et al [IceCube coll.], proceedings of the 31st ICRC, Lodz, 2009, arXiv:1004.2093 [astro- 
ph.HE]. 

[25] A. Achterberg et al. [IceCube coll.], ApP 26 (2006) 282. 

[26] J. Braun, M. Baker, J. Dumm et al, ApP 33 (2010) 175. 

[27] M. Baker, in these proceedings. 

[28] L. Demiroers, proceedings of Neutrino 2010, Athens, 2010. 

[29] M. Chernyakova, A. Neronov and M. Ribordy, proceedings of the 2009 Fermi Symposium, eConf 

C091122, arXiv:0912.3821 [astro-ph.HE]. 

[30] A. Neronov and M. Ribordy, PRD 79 (2009) 043013. 

[31] H. R. Christiansen, Phys. Scripta T127 (2006) 82. 

[32] G. E. Romero, H. R. Christiansen, M. Orellana, ApJ 632 (2005) 1093. 

[33] M. Orellana, G. E. Romero, Astrophys. Space Sci. 309 (2007) 333. 

[34] G. E. Romero, M. Orellana, A. T. Okazaki et al, Int. J. Mod. Phys. D17 (2008) 1875. 

[35] D. F. Torres, F. Halzen, ApP 27 (2007) 500. 

[36] J. Albert et al, Science 312 (2006) 1771; H. Anderhub et al, ApJ 706 (2009) L27. 

[37] V. A. Acciari et al, ApJ 679 (2008) 1427. 

[38] T. Tanaka, Y. Fukazawa, T. Mizuno et al, AIP Conf. Proc. 921 (2007) 391. 

[39] D. F. Torres, S. Zhang, J. Li et al, [arXiv: 1007.2272 [astro-ph.HE]]. 

[40] A. Smith, P. Kaaret, J. Holder et al, ApJ 693 (2009) 1621. 

[41] P. C. Gregory, ApJ 575 (2002) 427. 

[42] J. Casarcs, I. Ribas, J. M. Paredes et al, MNRAS 360 (2005) 1105. 

[43] V. Zabalza, J. M. Paredes, V. Bosch-Ramon, [arXiv:1011.4489 [astro-ph.HE]]. 

[44] J. M. Paredes, J. Marti, M. Peracaula and M. Ribo. A&A 320 (1997) L25. 

[45] A. A. Abdo et al, ApJ 701 (2009) L123. 

[46] R. Abbasi et al. [IceCube coll.], PRD 79 (2009) 062001; http://www.icecube.wisc.edu/science/data. 

[47] J. Abraham et al. [Pierre Auger coll.], ApP 29 (2008) 188. [Erratum-ibid. 30 (2008) 45]. 

[48] P. Abreu et al. [Pierre Auger coll.], ApP 34 (2010) 314. 

[49] V. Berezinsky and G. Zatsepin, Phys. Lett. B 28 (1969) 423. 

[50] R. Engel, D. Seckel & T. Stanev, PRD 64 (2001) 093010. 

[51] M. Ahlers, L. A. Anchordoqui, M. C. Gonzalez-Garcia et al, ApP 34 (2010) 106. 

[52] L. A. Anchordoqui et al, PRD 76 (2007) 123008. 

[53] R. Abbasi et al [IccCubc coll.], PRD 82 (2010) 072003. 

[54] K. H. Kampert [Pierre Auger coll.], AIP Conf. Proc. 1085 (2009) 30. 

[55] M. Ackermann et al, J. Geophys. Res. Ill, D13203 (2006). 

[56] G. Askaryan, JETP 14 (1962) 441; G. Askaryan, JETP 21 (1965) 658. 

[57] G. Askaryan, Sov. J. Atom. Energy 3 (1957) 921. 

[58] G.A. Askariyan, B.A. Dolgoshein et al, NIM 164 (1979) 267. 

[59] J.G. Learned, Phys. Rev. 19 (1979) 3293. 

[60] S. Bevan et al, ApP 28 (2007) 366. 

[61] L. G. Dedenko et al, Bull. Russ. Acad. Sci. Phys. 61 (1997) 469. 

[62] E. Zas, F. Halzen & T. Stancv, PRD 45 (1992) 362. 

[63] I. Kravchcnko et al, ApP 20 (2003) 195; I. Kravchenko et al, PRD 73 (2006) 082002; D. Z. Besson, 

J. Phys. Conf. Scr. 81 (2007) 012008. 

[64] P. Chen and K. Hoffman, preprint arXiv:0902.3288. 

[65] S. Boeser et al, J. Phys. Conf. Scr. 81 (2007) 012016. 

[66] R. Abbasi et al, arXiv: 1004. 1694 [astro-ph.IM]. 

[67] B. Price, J. Geophys. Res. Ill (2006) B02201. 

[68] T. Meures, L. Paul and M. Ribordy, proceedings of the 4th ARENA conf., 2010, Nantes, 2010, 
arXiv: 1010.2313 [astro-ph.IM]; M. Podgorski and M. Ribordy, proceedings of the 31st ICRC, Lodz, 
2009, arXiv:1001.3963 [astro-ph.IM]. 



10 



[69] R. Abbasi et al. [IceCube coll.], ApJ 718 (2010) L194. 
[70] P. Desiati, f. t. IceCube coll., arXiv: 1007.2621 [astro-ph.IM] . 

[71] K. Nagashima, K. Fujumoto and R. M. Jacklyn, J. Geophys. Res. 103 (1998) 17429; D. L. Hall et al, 
J. Geophys. Res. 104 (1999) 6737; A. Amcnomori et al, Science 314 (2006) 439; J. L. Zhang et al, 
proceedings of the 31st ICRC, Lodz, 2009; G. Guillian et al, PRD 75 (2007) 062003. 

[72] A. A. Abdo et al, ApJ 698 (2009) 2121. 

[73] A. H. Compton and I. A. Getting, Phys. Rev. 47 (1935) 817. 

[74] A. Lazarian, P. Desiati, ApJ 722 (2010) 188. 

[75] R. Abbasi et al. [IceCube coll.], ApJ 710 (2010) 346. 

[76] E. Waxman and J.N. Bahcall, PRL 78 (1997) 2292, 

[77] J. Ahrcns et al. [AMANDA coll.], ApP 16 (2002) 345. 

[78] F. Halzcn, J. E. Jacobscn, E. Zas, PRD 53 (1996) 7359. 



